Nonrelativistic Particle in Free Random Gauge Background 



Michael Engelhardt* 
Department of Condensed Matter Physics 
Weizmann Institute of Science 
Rehovot 76100, Israel 



The problem of a nonrelativistic particle with an internal color degree of freedom, with and 
without spin, moving in a free random gauge background is discussed. Freeness is a concept de- 
veloped recently in the mathematical literature connected with noncommuting random variables. 
In the context of large-N hermitian matrices, it means that the the multi-matrix model considered 
contains no bias with respect to the relative orientations of the matrices. In such a gauge back- 
ground, the spectrum of a colored particle can be solved for analytically. In three dimensions, near 
- - - zero momentum, the energy distribution for the spinless particle displays a gap, while the energy 

, distribution for the particle with spin does not. 
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I. INTRODUCTION 



Much of the complexity of strong interaction dynamics derives from the non-commutativity of the gauge fields, 
encoded in their matrix character. It leads to the gluonic self-interactions, which generate phenomena such as 
asymptotic freedom, and presumably also confinement due to the constriction of chromoelectric fields into fiux tubes. 
Matrix-valued degrees of freedom are notoriously difficult to handle, even in cases where one has managed to eliminate 
I the space-time dependence of the matrices. One such example may be the large- A^c limit of QCD {Nc denoting the 
. number of colors) |Q] Q . There one argues that gauge field path integrals should be dominated by a large- A^c saddle 
' point, and that there is a gauge where this saddle point, commonly referred to as the "master field", is constant in 
^ : space and time. This is plausible in view of the translational invariance in space-time of the underlying action and 
measure. 

, In the case of large- A^c" QCD, it is not even known how to formulate a potential purely in terms of space-time 
■ independent matrices which will reproduce the master field; however, also if one posits some potential of one's choice, 
O i] such simplified matrix models generically become intractable as soon as more than one matrix degree of freedom 
1^ I is involved, i.e. as soon as non-commutativity is allowed to play a role. There are notable exceptions, such as the 
1^ ' Itzykson-Zuber integral [||, which lies at the heart of the Kazakov-Migdal and Penner models [Q 0. In these models, 
J> , due to the special type of potential involved, some selected observables can be calculated. 

Recently, there has been some advance in the mathematical literature concerning a type of non-commuting variables 
called free random variables it has furthermore been shown (loc.cit.) that in the large-N limit, independent 
' hermitian matrix models describe variables of this type. These are models governed by a partition function of the 
form 

Z ^ [l[[dA,]cxpUNj2T^ViAM (1) 

i=i \ i=i / 

where the integration is to be carried out over the real and imaginary parts of the matrix elements of the Ai, subject to 
the constraints of hermiticity. The potential, consisting of a sum of terms each only involving one of the matrices, thus 
only determines the eigenvalue distributions of the Ai but there is no bias with respect to the relative "orientations" of 
the matrix degrees of freedom; each orientation is equally probable. This restriction is quite strong; most interesting 
matrix models, presumably including the one describing the master field of QCD, contain interactions between the 
different matrices. It should be noted, though, that the constructions used in the theory of free random variables may 
be extended to more general models and have recently led to some increased understanding of the functional analytic 
spaces on which the master field may be formally defined § [§ [|. 



* email: hengel® wicc . weizmann. ac .il 



1 



However, even noninteracting matrix models constitute interesting laboratories for the study of phenomena asso- 
ciated with non-commutativity. Initially, it is well known how to derive the eigenvalue distributions of the involved 
matrices from their respective potentials. The eigenvalue distribution is related by a Hilbert transform to the deriva- 
tive of the potential; i.e. if the hermitian matrix A is governed by a potential term TtV{A), then the eigenvalue 
distribution (normalized to unity) is determined by 

lv'iz)= I d\p{X)^ (2) 

Then, however, as soon as one attempts to calculate eigenvalue distributions of objects composite in the different 
matrices contained in the model, one is faced with their non-commutative nature. Here, the methods developed in 
exhibit their full power; they enable one to convolute the distributions of different free random variables in a quite 
general and systematic manner, not restricted to a very special form of the objects considered. 

In this paper, two not altogether trivial cases of problems which can be solved completely analytically with the 
techniques developed in are exhibited. These are the problems of nonrelativistic particles with an internal color 
degree of freedom, with or without spin, moving in a free Gaussian random gauge background. The objective of 
this work is to develop some intuition regarding the phenomenological consequences of freeness in the context of 
large-N matrix degrees of freedom. This can be viewed as complementary to the recent more formal investigations 
mentioned above concerned with the formal definition of master fields and the spaces on which they operate. Here, 
by contrast, the emphasis lies with the study of some exactly solvable model systems which may serve as paradigms 
for more complicated realistic cases. The examples treated here, while sufficiently simple to be of pedagogical value, 
nevertheless seem generic enough to be relevant for applications. Some physics questions for which they may be of 
relevance will be mentioned at the end of the paper. 

It is to be expected that developing some insight into the phenomena resulting from freeness also is of some 
value from the point of view of the quest to understand interacting matrix models beyond the few existing solvable 
cases. With the advent of techniques allowing one to deal with free variables, one controls the two limiting cases of 
vanishing and infinitely strong coupling between matrix degrees of freedom: A system of matrices whose orientations 
are completely aligned due to some strong coupling is as simple as a one-matrix model, since all matrices can be 
simultaneously diagonalised; on the other hand, a system of matrices whose relative orientations are subject to no 
bias at all can now be treated as well. Knowing both limits should provide valuable qualitative predictions for the 
behaviour of models even with nontrivial couplings. 



II. PARTICLE WITHOUT SPIN 

Consider the Hamiltonian 

D 

Ho = ^(A:, - A,)^ (3) 

1=1 

This is the Hamiltonian of a spinless nonrelativistic particle with color in D dimensions, minimally coupled to a 
background vector field A. Each of the components of A is taken to be an x A^ hermitian matrix (where A^ oo), 
constant in space and time. Due to this last condition, the Hamiltonian is diagonal in momentum space and the 
solutions can be classified according to the momentum components ki. For given k and A, there are A^ different 
modes of propagation, i.e. color eigenvectors, for the particle. Now, consider the vector field not as fixed, but to be 
averaged over an ensemble described by an independent matrix model. In other words, any observable O should be 
averaged as 

= i / Yl[dA.]0{A)exp (-A^f]TrF(A,) ) (4) 
•' j=i \ j=i / 

Note that now the model is invariant under arbitrary unitary rotations of the variables Ai. 

In short, the problem proposed is the following: Given the eigenvalue distributions of the Ai ( which can be easily 
found from the potential V{Ai), as mentioned in the introduction), find the eigenvalue distribution of Hq for arbitrary 
fcj. Since the model contains no bias with respect to the relative orientations of the Ai, the Ai constitute a family 
of free random variables, and the problem is a fairly straightforward application of the techniques developed in [pj. 
This will serve as a warmup for the more involved case of the spinning particle, where one will have to deal with more 
exotic situations, such as the free random variables themselves appearing as components of matrices. 
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To treat a concrete model, the potential V{Ai) must be specified. If one wishes to preserve invariance of the model 
under spatial rotations without destroying the freeness, one is forced to choose a Gaussian potential, V{A) = (3A^. 
The eigenvalue distribution for the matrices Ai induced by this potential is the well-known semicircular distribution 

Pa{\) = -^Va2-A2 (5) 
7ra 

where = 2//3. 

As a first step in finding the eigenvalue distribution of the Hamiltonian, one easily finds the eigenvalue distribution 
of each of the summands {ki — Ai)"^ of Hq: If one has pa(A)c?A eigenvalues of Ai in the interval [A, A + dX], then the 
number of eigenvalues of {ki — AiY to be found in the interval [(fc — A)^, (fc— A)^ — 2(fci — A)ciA] is {pa{X)+ pa{—^)) dX. 
Substituting e = {h — A)^, one thus finds p(^ki-Aiy {^)de eigenvalues of {ki — Ai)^ in the interval [e, e + de], where 

P{ki-Ai)^{e) = -^{PAiV^ + h) + pa{V^ - ki)) (6) 

Of course now there are only positive eigenvalues, e > 0. 

It remains to convolute the eigenvalue distributions of the different summands {ki — Ai)'^ . These are hermitian 
random variables, with eigenvalue distributions given by (^), and freely rotating with respect to one another, i.e. 
there is no bias in the relative orientations in color space of the variables when taking the ensemble average over the 
Ai. To convolute these eigenvalue distributions, one must compute their so-called R-transform The R-transform 
plays the same role for adding noncommuting free random variables as the logarithm of the Fourier transform docs for 
adding ordinary commuting random variables: Convolution becomes ordinary addition, i.e. adding the R-transforms 
of the eigenvalue distributions of the summands gives the R-transform of the eigenvalue distribution of the sum. 
The R-transform is defined by the following algorithm: First, find the generating function of the moments of the 
distribution: 

oo 

/oo 
dAp(A)A" (8) 
-OO 

which can be analyticaUy contmucd to the entire upper complex half-plane by summing the geometric series: 

G(0 = /^rfA|^ (9) 
This continuation is important for the inverse problem of going from the generating function back to the distribution 

p(A) = limImG(A + i77) (10) 

Next, find the function inverse to G, i.e. K such that K{G{C)) = £.■ Then, the R-transform is given by R{z) — 
K{z)-l/z. 

It turns out that this program can still be carried out analytically for the program at hand; one is led to algebraic 
equations of at most fourth order and thus just remains within the realm of solvability. The generating function 
corresponding to the distribution (||) is easily calculated to be 

GdO = f"^d,_L_PA(^^ + k^)+jA{V-^~h) 

de^^^^ (12) 







^In the following, a factor N will be taken out of eigenvalue distributions, i.e. pA is normalized to unity. 
^This problem has recently also been treated from a more physical point of view in 



3 



de - — ■= 

-1 {V^+k^ 



h + ae) 



(13) 
(14) 



Now, to find the inverse Ki of Gi, one must solve (|lj) for ^ in terms of Gi. By always bringing the terms containing 
no square roots to the left hand side and squaring the resulting equation (this procedure must be applied twice) , one 
arrives at an equation with no square roots; in this equation, quadratic in ^, the constant term vanishes, the trivial 
solution ^ = is spurious, and thus one arrives at 



R,{z) + - = K,{z) = ^{z) ^ - + 



4fc? 



a^z (2-a2z)2 



(15) 



The R-transform of the eigenvalue distribution of the Hamiltonian is now simply the sum of the individual R-transforms 
corresponding to each spatial component: 



D 



4fc2 



j=l 



A-a'^z (2-a2z)2 



(16) 



where fc^, the squared Euclidean length of the vector k, has been introduced. Note that the h dependence of 
the individual R-transforms Ri could have already been predicted from the additivity of the R-transform and the 
rotational invariance of the Hamiltonian^. From J^iLi fi^i) — -^(1^1) it already follows that / must be of the form 
f{ki) = c + dkf, as is manifest in (p^). 

It remains to go from the R-transform Rhq back to the eigenvalue distribution of the Hamiltonian. This is the 
more arduous task. Inverting Khq{z) = Rh„{z) + 1/z corresponds to solving the equation 



4(2 - a^z)^ + {D- 1)0^ z{2 - a^zf + Ak^z{A - a^z) = Cz(4 - a'^z){2 - a^zf 
for z in terms of ^. This is a fourth order equation for z, which via the substitution 



D-1 



can be brought to the normal form 



-f- (3iV - 6Af2)a;2 -|- (SM^ - mM)x + 3NM'^ - 3M^ + L/4 = 



where the abbreviations 



L = 



64fc2 



M 



D - 1 



N 



3a2^ 



1 



D - 1 



n2 



4 
3a4 



(17) 



(18) 



(19) 



(20) 



have been introduced. The analysis of such an equation is quite involved, and for this reason is relegated to the 
appendix. The final result for the eigenvalue distribution of the Hamiltonian (||), which according to (^0|) is essentially 
the imaginary part of the solution of (p^ , can be stated as follows: 



PHo iO 



1 

2^ 



0{~d){i ~ e{~b2)0{h)) 



where 



2/2 



2Re (we2"/3) , 
P ' 



2N + 4M2 

i\l ~d 
-3N^ 



6N 



L 
12M^ 



+ (1 - 20(4M3 - 3iVM))/^ 

ys = 2Re {ue'^'^'^^) ~ 2N + AM"^ 
^=27 + T 

q = -2^^'"^ + 2NL - AAPL 



L 



(21) 

(22) 
(23) 

(24) 
(25) 



^Remember that the original Gaussian potential governing the distribution of the vector field components is invariant under 
rotation of the vector field in the spatial indices. 
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The solution p^/7jg-(A) = 2ApiJo(A^) is plotted (this representation avoids the trivial square root singularity at zero 
eigenvalue, especially further below in the case of the spinning particle) for various values of the free parameters k/a 
and D in figures (0)-(|). Note that when comparing different Z?, there are two distinct sensible ways of choosing 
the units. On the one hand, one may take the width a of the distributions of the vector fields as the unit of the 
square root of energy, as one does when comparing different k (cf. figures (Q) and (^). However, one may also rescale 
a? = a? / D and use a as the unit of the square root of energy. This corresponds to keeping the norm of the vector 
field, i.e. Tr Af, constant as one varies D (just as one automatically does in the case of the momentum k by 

only talking about its modulus), cf. figure (||). 

To characterize the behaviour of the solution more precisely, it is useful to evaluate various moments and limits. The 
first few moments of the distribution pn^ can be obtained without knowledge of the full solution (^l|). It is sufficient 
to expand the R-transform ( |l6| ) in a power series in z and invert the corresponding scries Kho{z) = Rhq{z) + l/z to 
the desired order, yielding the first few terms of G(^) expanded in powers of 1/^. The moments are then by definition 
the coefficients appearing in G'(^), yielding 

>^o)^^ + fc^ = 5 + fc^ (26) 
— Tr Hi) - ( — Tr Hq ) ^ + a^k^ ^ + 27 

With the help of this information, one can obtain the behaviour of the distribution as /c — > 00 or 13 — > 00. 

Consider large D, and use as the unit of energy. The distribution is centered around a^D/A, with a width 
behaving as \/Z3. Substituting therefore 

C=^+7V^ (28) 

in ( |2l|) , and keeping only leading pieces in I/a/D, one obtains that 62 < always, and 61 < for 7 e 
[-oV#T^/2,aVFT^/2]. Also, 

218(472 _a4)p 1 , ^ 

Lao JT. (29) 



i.e. d < for 7 G [—0^/2, 0^/2]. The condition on d is more stringent than the one on bi and thus defines the support 
of the distribution. Furthermore, one has 

M2 - 37V/4 = - 4 (a^ (a^ - fc^ ) - 47^ ) — (30) 
a" D 

lyal = ^ (a2(a2 + k^) - 47^ + 2aky/a^ - 472) 1 (31) 



12/31 



^ ( a\a^ + fc2) - 472 - 2ak^a'i - AjA - (32) 



lysl = Vly2l + ly3|±2Vly2||y3| (33) 

: ^y/a^a^ + fc2) - 472 ± |a2(a2 _ ^2) _ 472]^ (34) 



Pff„(e) = — ^V«*-472 (35) 



and therefore 



In other words, for D — s- 00 with fixed k and a one obtains a semicircular distribution centered at a^D/A and with 
radius 

On the other hand, if one chooses to use a'^ as the unit of energy, the limiting behaviour as D — > 00 can already be 
read off from the first two moments. The width of the distribution vanishes and one obtains a 5-peak at a2/4 + fc2. 

Finally, taking the limit /c — s- 00 can be done in complete analogy to the limit D — > 00 above. However, it is not 
necessary to explicitely go through this: For fc ^ 00, one can approximate the Hamiltonian by Hq ~ + 2 X^iLi kiAi. 
This means that one simply has to additively convolute the original semicircle distributions of the Ai , with the radii 
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scaled by a factor 2ki , respectively. It is well known Q that the semicircle distributions are closed under additive free 
convolution and that the radii add like components of an Euclidean vector. Therefore one immediately obtains the 
result that, for fc ^ oo, the eigenvalue distribution of Hq approaches a semicircle centered at and with radius 2ka. 

The physical interpretation of this result is clear: For k oo, the quantum mechanical propagation of the particle 
is dominated by the semiclassical straight trajectories in the direction of k. The particle does not see the transverse 
directions and behaves just as it would in a one-dimensional world; correspondingly, the dimension D does not enter 
the limiting semicircular eigenvalue distribution as fc — > oo. 

On the other hand, when k is finite, the particle does explore the dimensions available to it to some extent. For 
large D, the width behaves as a^^/D — o^/^/D (and thus does not tend to a finite value when one scales the vector 
field to have constant norm as D is varied). The semicircular nature of the distribution in the limit of large D is 
understandable as a manifestation of the central limit theorem: For a large number of variables (matrices) , the sum of 
the random matrices (after separating off the first moments) behaves as if it was distributed according to a Gaussian 
weight, i.e. as a semicircle 0. 

Most interesting is the soft region around k — 0. For fc = 0, the solution (^) simplifies considerably. In this case, 
one has 

d = (36) 
62 = 6iV - 12M^ (37) 
&i = {3N - 4:M^){3N - 12M^) (38) 

and therefore only a nonvanishing eigenvalue distribution when N > 4M^/3. In this case, one has ■y/|y2| — VTysJ — 
V 3N — 4Af ^ and therefore 

p^^ (^) = -Q(ZN- 4M2) ^JzN - 4M2 (39) 

TT 

The edges of the distribution are determined by 



= 3A'-4„» = -(£-i£ + Hffi±il-l (40, 

44 a c, " 



leading to ^ = a^^sfD ± 1)^/4, i.e. the support of the distribution is the interval {a^^s/D - l)"^ /A, a'^{\^ + 1)^/4]. 
Thus, in two or more dimensions, the eigenvalue distribution has a gap, even for k 0, as already evidenced in 
figures (0)-(^. Of course, for _D = 1 one just has for Ha the distribution (^), which behaves proportionally to l/VJ 
as ^ ^ 0. Note also that the above result on the support of the distribution implies that, at fc = 0, the spectrum of 
the square root of Hq always has a support of length a, for any dimension. 

Formally, the emergence of an energy gap does not seem particularly surprising in view of the fact that the 
Hamiltonian (^ is a sum of squares. Consider the more familiar case of ordinary (commuting) random variables, 
taking positive values, and with a finite probability density when one approaches zero. Typically, the distribution 
of a sum of two such variables will vanish linearly at zero, because sampling the two variables, it is improbable to 
simultaneously find two very small values. There is actually a matrix realisation of this, namely two random matrices 
with the aforementioned type of distributions for the eigenvalues, but not freely rotating over the whole U{N) group 
with respect to one another; instead, constrained to commute, i.e. only all permutations of the eigenvalues are 
allowed. Presumably to realise this, one would have to introduce a strong interaction term between the two matrices, 
proportional to the square of their commutator. 

In the present case, by contrast, one is allowing a much larger class of relative orientations between the matrices, 
namely an arbitrary U{N) rotation. This larger invariance class evidently enhances the effect of making the occurrence 
of a small eigenvalue improbable in the sum of the matrices. For free random variables, the effect is so strong that the 
resulting distribution not only vanishes as one approaches zero, but vanishes on a whole interval between zero and a 
finite value. 



III. PARTICLE WITH SPIN 



The Hamiltonian of a particle with spin in the same gauge background as was considered in the case of the spinless 
particle is the 2N x 2N matrix 

Hs = {a{k - A))^ = {k - Af + ia{A X A) (41) 
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where the ai denote the Pauh matrices; the number of space dimensions has now been speciahzed to D = 3. In the 
last term on the right hand side, one recognizes the couphng of the spin to the nonabehan magnetic field; the latter 
contains no derivative term since the vector potential is taken to be spatially constant. 

In practice, it is not helpful here to separate the Hamiltonian into the spinless part and the spin-magnctic-field 
coupling. The crucial question when applying the free random variable techniques is whether one can cast the 
calculation into successive multiplications and additions of mutually free variables. E.g., given two free random 
variables X and Y, it is in general not straightforward to calculate the eigenvalue distribution of + XY; one 
may accomplish the multiplicative convolution in the second summand (how this is done is explained below), but 
the two summands do not rotate independently of one another and therefore the free convolution techniques do not 
apply. On the other hand, the eigenvalue distribution e.g. of X + XY is straightforward to obtain once one writes 
X + XY = X{l + Y). Thus, in the case of the Hamiltonian having solved the spinless problem is of no assistance; 
one must restart and consider the eigenvalue distribution of the object 

P^a{k^A) (42) 

(obtaining the distribution of Hs — at the end is trivial). 

As in the case of the spinless particle, the eigenvalue distributions of the Ai are taken to be semicircular with radius 
a; the harmonic potential generating these distributions is invariant under spatial rotations. Therefore, without loss 
of generality, one may take k in the 3-direction, k = (0, 0, k). Being able to do this will be crucial for the developments 
to follow. 

One is faced here with a new situation: Free random variables appearing as components of a matrix. Considering 
for the moment matrices of finite rank, the eigenvalues A of P for a specific realisation of the gauge matrices Ai are 
determined by the equation 

= dct(P- A) = dctgA (43) 

with the N X N matrix 

Qx=A3-k-X+{Ai+ iA2){A3 - fc + X)-\Ai - iAa) (44) 

Note that Qx has only half the rank of P — A; however, since Q\ is nonlinear in A, it can have the same number of 
singular points in A as P — A, as required. The second equality in ( ^ ) is only valid as long as A3 — fc + A has no exact 
zero eigenvalues. However, this happens only on a set of A of zero measure; it suffices here to consider the generic 
case where the above equality is valid. Also, the same consideration for P = {10-2) P{— 1(^2) shows that the ensemble 
averaged eigenvalue distribution of P is symmetric about zero if the eigenvalue distributions of the Ai are. Therefore, 
it will be enough to consider positive A. 

When one is faced with matrix equations nonlinear in a parameter such as A in the present case, it is customary to 
go to an equation in larger matrices linear in the parameter |]l3| . This is also familiar in the context of "supersymmetry 
tricks" . Here, it is not helpful; it just leads from Q\ back to the original matrix P. Instead, here the extended problem 

det(QA - Ai) - (45) 

will be considered, and at the end specialized to ^ = 0. In the large-iV, ensemble averaged language, the eigenvalue 
distribution of Qx, parametrically depending on A, will be obtained using free convolution techniques. This yields a 
distribution px (fJ-) which describes how many solutions of ( p5| ) one finds at a point in the X — fj, plane if one counts 
them as they occur on an increment dfi in /i-direction. In the end, this will have to be translated into how many 
solutions one finds if one counts them as they occur on an increment dX in A-direction; then, one may set fj, — 0. 

One further issue must be settled before the free convolution techniques can be applied: In general, the two 
summands Qi and Q2 appearing in Qx = Qi + Q2, 

Qi=A3-k-X (46) 
Q2 = (Ai + tA2){As -k + X)-\Ai - iA2) (47) 

do not rotate independently of one another, since they both contain ^3. This would thwart attempts to use free 
convolution. The special case considered here, however, can be rescued: The matrix A\ + 1^2, with A\ and A2 
hermitian, represents a general complex matrix, which can be reparametrized using the polar decomposition p4| 

Ai + iA2 = V^BVU'' (48) 
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where U and V are unitary and B is diagonal and positive. Integrating over all hermitian Ai and A2 corresponds to 
integrating over U and V with Haar measure and over the eigenvalues rj contained in _B, 



[dAi] [dA2] = J{rf)d^ri [dU] [dV] (49) 

where the Jacobian is 

N 

J{ri) = X{ri.X\{vf-ri^,f (50) 

i—l i<j 

To be precise, the integration over V is only over the quotient group U{N) /{U{1)^); this is no consequence, as will 
be explained below. Also, note that in (^8|), the product VU^ is usually written as a single matrix; however, since U 
is being integrated over with Haar measure, one can always pull out a factor V . 
In the new variables, the potential governing Ai and A2 becomes 

Tr V { Ai, A2) = PTt{AI+ Al) ^ f3TTB^ (51) 

and one has 

Q2 = V''BVUHA3~ k + X)-^UV''BV (52) 

Cast in this form, it becomes clear that the calculation of the Q\ eigenvalue distribution can be carried out using 
free convolution: The combination [/^(As — k + X)~^U rotates independently of Qi = A^ — k — X, and therefore the 
evaluation of the eigenvalue distribution of Qx indeed becomes a sequence of free convolutions if one additionally uses 
that the eigenvalue distribution of Q2 is the same as the one of 

Q2 = {V^BV)-'^Q2{V''BV) = U\A3 -k + Xy^UV^'B'^V (53) 

(up to regions of zero measure in B) . It is also clear that the restricted integration domain of V is of no consequence: 
For two variables to be free with respect to one another, it suffices to rotate one of the variables. The problem is in 
fact invariant under the additional rotations by V and the 1^-integration gives only a normalization factor. 

Finally, before applying multiplicative free convolution to obtain the eigenvalue distribution of Q2, it remains to 
solve the saddle point equation for the eigenvalue distribution of B^ , which is controlled by the Jacobian (|5^) and the 
potential (pl|). This standard procedure is carried out in the Appendix, yielding 



1 / — n 

PB2(77) = — J 0<r,<2a2 (54) 

Note that B itself is therefore distributed according to a quarter-circle of radius y/2a. 

As explained in connection with the spinless particle, additive free convolution of eigenvalue distributions is accom- 
plished by going to the corresponding R-transforms, which are additive; similarly, for multiplicative free convolution 
one defines the S-transform which behaves multiplicatively. The S-transform is defined as follows: Find again the 
moment generating function G{z) as in (^; solve for the function inverse to G{l/z)/z — 1, i.e. T(y) such that 

\G(-^)-l = y (55) 



Then the S-transform is given by 



T{y) \T{y) 



S{y) = i^T(y) (56) 



With the help of these techniques, one can now carry out the necessary convolutions leading to the eigenvalue 
distribution of Q\. Starting with the eigenvalue distributions (jsj) and (cf. (^),(|1)) 

P(A3_fc+A)-i iv) = -^PA Q + fc - (57) 
one obtains the corresponding generating functions 
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GiA,-t+^)-{z)= f dx-^i^— (59) 

J ~oo Z x-k+X 

/I 1 \ 

(60) 



dx\J cP- — I 1 H T-\ 

V zx-k+X-^J 



^ + k-X- ^{-^+k-Xf -a-^j (61) 

In both cases, T is determined by a quadratic equation after one has apphed a procedure analogous to the one leading 
to eq. (^5|). In the case of Tb^, the correct solution is easily picked out because there is a spurious zero solution; in 
the case of T(^a^_}^^x^-i ^ the correct solution can be picked out by considering the limiting behaviour as a — > and 
using that the S-transform of the unit matrix multiplied by a constant w is 1/w. Thus one arrives at the S-transforms 

Ss^iy)^^-^ (62) 



v3-...)-^(.)-^ i+ji-^?;^ (63) 



{X-kf 



(64) 



From this, one obtains that the function 

^^-l.T-(.).. (65) 

satisfies the equation 

z^a'*(l + uf - 2za'^{X - k)u{l + u) + a^u^ = (66) 

and consequently, Gq^ [z) = v satisfies 

a^v^z - 2a^{X - k)v{vz - 1) + a^ivz - 1)^ = (67) 

There is no need to solve for Gq^ at this point; the eigenvalue distribution of Q2 must still be additively convoluted 
with the one of Qi. Therefore, one can solve ( |67|) directly for z = Gq^{v), 



- = i - + A - A- (A - - ^) - ^ (68) 

where the relevant solution has been picked out by using, in the limit a — )■ 0, that Q2 ^ and therefore the 
R-transform Rq2{v) = Kq^iv) — 1/v — )■ 0. Combining this with the R-transform corresponding to Qi (cf. (46)), 



RQ^{^) ^ - k - X (69) 



one finally obtains 

1 a" 



V 



G-\[v)^KQM = ---r^'^k-{X-k)J[l--——\ (70) 



4 ' V 2{X-k)J [X-kf 



This determines Gq^ and thus ultimately p\{fJ.) ~ —lmGQ^{fi)/'K. However, before embarking on the arduous task of 
inverting ( |70| ) , it is advantageous at this stage to consider how p\ (fi) in the end determines the eigenvalue distribution 
pp{X) of P (cf. (|4^)). Considering for the moment finite matrices, the solutions of equation ( |45| ) for a specific 
realization of the gauge matrices Ai define continuous trajectories in the A — ^ plane, cf. figure (0). Asymptotically, 
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there are trajectories in the vicinity of the hnes defined hy X + k = and A — fc = 0, where "vicinity" means at a 
finite distance of roughly up to the semicircle radius a. The quantity one is ultimately interested in is the number of 
trajectories one crosses if one marches from A = cx) in along the A-axis to some point Aq; this is essentially the integral 
over pp{X). Consider the conjecture that this is the same as the number of trajectories one crosses marching from 
(/i, A) — (cx3,cx3) in A-direction to (/i, A) = (oo,Ao), and then in /i-direction to (/i, A) = (0, Aq), cf. figure There 
is one point which must be clarified before this statement can be accepted as true: In general, it might happen that 
a trajectory intersects the integration paths described above more than once, and then the two paths may count a 
different number of trajectory crossings, as displayed in figure in the right-hand graph. It will now be argued that 
this is not possible. To begin with, note that 

gA(fc)-M = QA+M/2(fc + Ai/2) (71) 

This means that solving the extended problem (^) with general fj, is in fact equivalent to solving the original problem 
( ^ ) with shifted k and A (conversely, k could be absorbed into A and ^ ) . Now, the solutions A of the original problem 
det(P — A) = are continuous in k, which can be seen as follows: Since the determinant of a Hermitian matrix is real, 
the characteristic polynomial det(P— A) is real-valued for real A; its coefficients are real polynomials in k. Therefore, 
the graph of det(P — A) varies continuously with k. Now, if one wanted to "annihilate" two zeros of det(P — A) as 
k is varied, say by having a minimum of det(P — A) cross from below the A-axis to above the A-axis, on would have 
to "create" two other zeros somewhere else on the A-axis in a similar fashion, because det(P — A) must always have 
2N zeros {2N denoting the rank of P). However, at the point in k where this discontinuity takes place, one would 
thus have more than 2^^^ zeros, which is impossible. Therefore, the eigenvalues are indeed continuous as k is varied. 
According to (|7l|), this immediately also implies that in the extended problem det(Q\ — /x) = 0, as /i is varied, the 
trajectories X{fJ.) must be continuous. 

Consider now the possibility that a trajectory in the X — fi plane crosses a line of constant fi more than once, as 
in figure (^ in the right-hand graph. This would imply a discontinuous dependence A(/i) and can therefore not arise. 
In conclusion, this shows that marching along any line parallel to the A-axis in the A — /i plane, one can cross any 
eigenvalue trajectory at most once, as in figure (^ in the left-hand graph. 

One could argue in a similar way for paths in the /^-direction, except for having to be more careful due to the poles 
in Q\. There is an easier way to handle paths in the /x-direction once one has treated the case of the A-direction: Since 
one crosses each trajectory at most once when integrating in A-direction, it is necessary to cross all of them to saturate 
the normalisation condition in A-direction. Then, however, one must also cross all of them if one marches in straight 
lines between the points (/i. A) = (cxd,cxd) — + (cxd,Ao) — > (— oo,Ao) — > (— oo,— oo). In other words, in this way one 
counts at least 2N trajectories, and more if there are multiple intersections when marching in /i-direction. However, 
one can easily check at the end using the large- iV, ensemble averaged result that the integral is exactly 2 (remember 
that the factor N has been scaled out everywhere in the eigenvalue distributions) and therefore realisations of the 
Ai which generate multiple intersections when marching in /i-direction give vanishing contributions to the large-iV 
eigenvalue distributions. 

In essence, the basic property which makes these arguments work is that the trajectories in the X ~ fi plane can 
be interpreted in two different ways. Obviously, when considering the extended problem (^), one is considering the 
behaviour of the eigenvalues /i under changes of the parameter A. On the other hand, due to the equivalence (^l]), one 
can also interpret the trajectories as describing the behaviour of the eigenvalues A under changes of the parameter /i 
(up to the trivial additional linear shift in A). These two interpretations taken together ensure the monotonicity of 
the trajectories, i.e. that they can cross each line of constant A or /i only once. 

As a result of these observations one can now indeed relate p\{fJ.) to pp{X): 

/•oo /"OO /"OO 

/ dXpp{X)= / dppXo{p)+ / dXpextra{X) (72) 

which by differentiation with respect to Aq gives pp. The second contribution on the right hand side arises as follows: 
For |A — fc| < a, there is a nonvanishing density of trajectories as /i — > oo due to the poles in Qx (cf. figure (^). These 
trajectories asymptotically become parallel to the /i-axis, and therefore pxifJ-) vanishes (as it must due to normalisation) 
as /i — > oo, since px{p)dp measures the number of trajectories encountered when marching in /i-direction. On the 
other hand, marching in A-direction, one crosses these trajectories perpendicularly. For a concrete realization of the 
Ai, they are determined as follows: In the complete problem det((3i + Q2 — fJ-) = 0, one can neglect Qi for /i — > 00, 
i.e. one is looking for the infinite eigenvalues of Q2; these occur exactly when A3 — k + X has zero eigenvalue, i.e. 
when —A is an eigenvalue of A3 — k. Therefore Pextra{X) is just the (shifted) original semicircle distribution, 

Pextra{X) ^ PA{k - X) (73) 
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One can therefore now specify pp(A): 



PpW 



d_ 



d^pxip) + pAik~ A) 



1 d f°° 

-Im^y^ d^GQ^ip.) + PA{k - X) 



1 d '•0 

TT 0\ 



d 





d 


— Im 




TT 








— Im 


f 


TT 


JGq 



dxx—Gn^ix) + pA{k ~ A) 





1 d f 

-Im— dx GqI (x) +pA{k- A) 



-Im 



2(A-fc) J (A-fc)2 



(A-fc)W 1 



+ PA{k-X) 



(74) 
(75) 
(76) 

(77) 

(78) 

(79) 



(in this derivation it has been used several times that ImGQ;^(0) never diverges, which will be verified below). The 
last two terms in ( [79| ) cancel, and using = Gq^(Gq^(0)) together with (|70|), one can additionally eliminate the 
square root to arrive at 



PpW = 



1 



a^TT VGq,(0) 



-tGqAO) ---^ImGQ,(0) 



(ReGQ,(0))2 + (ImGQ,(0)) 



a 

2 4 



(80) 



It only remains thus to solve (70) for Gq^(O). For v = Gq^(O) in (70), the left hand side vanishes and the square 
root can be eliminated by squaring the equation. Eliminating the cubic term in the resulting quartic equation by the 
substitution 



one arrives at 



with 



X 4A 
^ = - + 



+ C2X^ + CiX + Co = 



C2 = 


-3< 




64 


Cl = 




^27 




16 


Co = 


^27 



{2L^ + 12KL^ + ISK^L + 3L- 9K) 

/QOJi^rS I ^Q(<^2r2 I Qr2 AQli^r I ( 



(81) 

(82) 

(83) 
(84) 
(85) 



given in terms of the dimensionless K = k/a and L — X/a. The analysis of this equation is relegated to the appendix; 
the final result for Gq^ (0) can be stated as follows: 



ImGgjO) 



2a 



0{dM^,~^)\ 



AL 



ReGQ^(O) = —{20(9K - 2L-^ - 12KL^ - 18K^L-3L) - l)Vyi+ — for d > 



3a 



(86) 
(87) 



where 



2/1 



62/3 


y2 = e" 


62/3 


62-- 


. 1/3 




i 


u' = - 



+ e^"/^!*' - 62/3 



^( 
3 ^ 

M 
3 



-1/3 



(89) 
(90) 
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and 



d=^ + ^ (91) 
27 4 ^ ' 

256 

p = — (2 - 3K^ - 9K^ - UKL - UK^L + + 2K^L^ + AKL^ - L^) (92) 
4096 

q = (7 - + 27 + bAK^ + A8KL + UAK^L + lOSK^L + 12L^ 

729 

+66K^L^ + ISK^L^ - 4:8KL^ - 56K^L^ + 3L^ - 6K^L^ + \2KL^ - 2L^) (93) 

The solution p^jj-^{\) = 2Ap//^,(A^) = 2pp{X)\xyo is plotted for various values of the free parameter K = k/a in figure 
(^. In the limit k —> oo, the coupling of the spin to the magnetic field in the Hamiltonian (41) can be neglected and 



the energy spectrum is therefore identical to the spinless case: Semicircular with radius 2ka, centered at k^. On the 
other hand, in the soft region around fc = 0, the two cases differ qualitatively: There is no gap in the spectrum at 
k = for the particle with spin, since the discriminant (|9l| ) is positive for fc = A = 0, giving a nonzero result for (|8^). 
Only at a certain nonzero value of k does the distribution detach from the origin. 

In order to characterize this spin-induced effect more precisely, it is useful to evaluate the first two moments of 
the energy distribution. Since in the present case, the full generating function Ghs corresponding to the distribution 
PHs is not explicitly available, one cannot directly read off the moments. However, they are easily calculated using 
the axioms of freeness Q, which the random matrices Ai obey. First, one trivially has that the first moments of the 
eigenvalue distributions of Hq and Hs are the same, 

\TrHs) = (i-TrHo)=k'+^-a' (94) 



2N / \N 

since the additional spin-magnetic field coupling in Hs obeys Tr {i(7{A x A)) = already in every fixed realization of 
the vector potential A. Then, one has 

— Tr ffA _ / J_Tr = / i_Tr H?\ - / — Tr Hr\ - / —Tr (n( A x A^al A x A\\\ 



.2iV^^-^^/-\2iV^'-^V ^\n^'''^/-\n^'''V -\^Tr(a(AxA)a(AxA))^ (95) 
(again having used the tracelessness of the cr-matrices) . Now, 

^Tr {a{A x A)o{A x A))^ = ^^Tr {A,A,A,A^ - A?A^)^ = -6 (^^ (96) 

Here, in the first equality it has been used that the matrices Ai obey identical distributions; in the second equality, 
the axioms of freeness for Ai and A2 have been utihzed along with ((l/A^)Tr A^) — and ((l/iV)Tr^f) = a? /A. 
Therefore, 

— TriJ|\ - / —Ti HsS = ( -Tr^o'^ - / -Tr^oV + -a^ = Pa" + — (97) 
27V ^/ \2N ^/ \iV V \^ / 8 16 ^ ' 

Thus, the effect of the additional interaction of the spin with the magnetic field is merely to broaden, but not to shift, 
the eigenvalue distribution of the Hamiltonian. It does this sufficiently strongly to make the energy gap arising in the 
case of a spinless particle disappear at small momenta k. Figure (^ compares^ P-JhZ P^/IT^I'^ '^'^^ different k. 

Evidently, it is possible for the particle to align its spin with the nonabelian magnetic field in ways which allow it 
to lower the energy associated with some of the modes of propagation to values near zero. There are sufficiently many 
such possibilities to make the gap in the energy distribution disappear for small momenta k. A calculation of more 
specific spin-color magnetic field correlation functions, which would give a more detailed description of this effect, will 
be foregone here. 



A factor two is divided out of Prjf^ in order to compare two distributions normalised to unity. 
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IV. SUMMARY 



In this work, the spectrum of a nonrclativistic particle with internal colour degree of freedom moving in a constant 
Gaussian free random gauge background was discussed. Both the case of a spinless particle and the case of a spin-1/2 
particle were considered. The limit of large momenta k is easily understood perturbatively; the energy distribution 
becomes semicircular. Also, in the spinless case for a large number of spatial dimensions, one obtains a semicircle as 
a consequence of the central limit theorem. On the other hand, the nonperturbative regime near k = displays very 
interesting features. The spectrum of the spinless particle exhibits a finite gap in two or more dimensions. Evidently, 
for free random variables, the suppression of small eigenvalues in sums of squares (such as the spinless Hamiltonian) is 
stronger than in the case of ordinary commuting random variables, for which only a linear vanishing of the distribution 
at zero occurs. By contrast, when the particle possesses a spin degree of freedom, the gap in the energy distribution 
disappears at small momenta. Formally, this is understandable since the Hamiltonian ceases to be a sum of squares; 
physically, the particle can align its spin with the background color magnetic field such as to lower the energies of 
some of the modes. This effect is sufficiently strong to obliterate the gap in the energy spectrum occurring in the 
spinless case. 

On a technical level, the treatment depended crucially on the fact that the Hamiltonians considered were block- 
diagonal in momentum space. Thus, it would be very hard e.g. to combine a gauge background such as the one 
used here with a spatial harmonic oscillator potential for the particle. On the other hand, one can easily envision 
introducing such a gauge background into a bag model. In such a framework, the gauge background would provide 
a mechanism for generating masses for the constituent quarks. It would be interesting to pinpoint the differences 
between a bag model with a random gauge background and a conventional bag model with constituent quark masses 
introduced by hand. 

A related issue is the question of chiral symmetry breaking in a free random gauge background. The chiral 
condensate can be related via the Casher-Banks formula (l5| to the spectrum of the Dirac operator [0 . The 
methods developed here for the 2x2 spin structure occurring in the Pauli Hamiltonian would seem to constitute a first 
step towards also evaluating the eigenvalue distribution of the 4x4 Dirac operator in a random gauge background. 

Finally, a different physical way to view the average over an ensemble of spatially constant gauge configurations is to 
identify the ensemble averaging with a domain or time averaging. If the gauge background through which the particle 
is propagating consists of domains of approximately constant color magnetic fields with random orientations, then the 
wave packet will feel different realisations of the gauge matrices as it evolves. Averaging over such a background can 
be replaced by averaging over one domain with all possible background gauge configurations if one argues the effects 
of the domain walls to be small. 



The author acknowledges useful discussions with R.Plesser, A.Schwimmer, and especially S.Levit. This work was 
supported by a MINERVA fellowship. 



APPENDIX A: SOLUTION OF FOURTH ORDER EQUATION DETERMINING EIGENVALUE 

DISTRIBUTION OF SPINLESS PARTICLE 



ACKNOWLEDGEMENTS 



The solutions of the equation 



X* + {3N - 6Af2)a;2 + (8A/^ - 6NM)x + iNAP - 3M^ + L/4 = 
are determined by the solutions of the corresponding cubic resolvent equation pj[ | 

+ b2y'^ + biy + bo = 



(Al) 



(A2) 



with 



62 = 67V - 12M^ 



bi = 9N^ - iSNM^ + iM'I* - L 



(8M^ - 6NM f 



(A3) 



as 



13 



1 



1 



1 



X2 

where in addition to the signs given, the signs of the square roots must be adjusted so as to guarantee 

8M^ - 6NM 



yivy2vy3 



(A4) 



(A5) 



In order to preserve continuity of the solutions in the parameters L, M, N, this latter condition in practice is handled 
as follows: Choose the square roots in the conventional way at som e initial point in parameter space where this is 
allowed. Points in parameter space where the right hand side of ( |A5D changes sign coincide with one of the solutions 
Hi of ( |A2|) touching zero; give the square root of this solution an additional minus sign when crossing such a point to 
maintain (A_5). Note that after two such sign changes one may have switched between solutions in (A4). 
The third order equation ([A2|) in turn can be simplified as follows. Via the substitution 



y 



it reduces to the normal form 



with 



and the associated discriminant 



2N + AM' 



+ pw + (? = 



- L 



-2N'^ + 2NL - UI^L 



^ " 27 ^ 4 " 4L2(m2 - 3N/4f + [AN^L + 2NL^){M^ - 3iV/4) 
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1/ 
27 



(A6) 



(A7) 



(A8) 



(A9) 



The latter form will be useful below. 

As a next step, it is necessary to pick out which of the solutions of the original fourth order equation is the relevant 
one. The only place where this question is easily answered is for ^ — > cxd. This will be chosen as the starting point 
and the solution will then be followed through parameter space using continuity. At ^ — s- cx), using the definitions 



(pO|), eq. (A2) simplifies to y{y — 4/a^)^ — and thus two of the solutions are yi^2 = A/a^. For the third solution, 
one needs to go to order 0(1/^^). Writing 7/3 = s/^^, and keeping only the leading order in 1/^ in (A2), one obtains 
s ~ {D — 1)^/4. Note thus that (A5) is fulfilled to leading order in 1/^ with the conventional choice of the square 
roots. Now, since the moment generating function GnaiC) must vanish as ^ — > 00, and 



D - 1 
4? 



(AlO) 



(cf. (18)), one easily picks out the solution X4 in (A4) as the correct one for ^ — > 00. 

Now, to follow the solution through parameter space, it is nece ssary to collect some facts about its behaviour. First, 
note that the constant term in the cubic resolvent equation (A2) is negative. By Vieta's theorem it thus follows that 
the product of its three solutions is always positive. Therefore, the solutions of the equation may behave in three 
qualitatively distinct ways: 

1. Three positive real solutions: d < A &2 < A 61 > 0. 

2. One positive real solution, two complex conjugate solutions: d > 0. 

3. One positive real solution, two negative real solutions: d < A (62 > V 61 < 0). 

where in the characterization of alternatives 1.) and 3.), Descartes' sign rule for the coefficients in ( [A^ ) has been 
used. 



Further to this, one can show that, for ^ > 0, changes of sign of square roots according to condition ( [A5) ca n only 
occur in regions of type 3.) in parameter space. The argument goes as follows: The right hand side of (A5) being 
zero implies 4M^ ~ 3N (note that M > for ^ > 0). In these circumstances, the cubic resolvent eq. (A2) simplifies 



to y — AM y — Ly = 0. Noting that i > for ^ > 0, this equation has one positive and one negative solution apart 



14 



from the one which is zero. One is therefore indeed in a region of type 3.) in parameter space; it is one of the two 
negative solutions which touches zero and whose square root therefore acquires a minus sign. 

Below, slightly more than this will be needed, namely that when the solutions go from a region of type 2.) to a 
region of type 3.) and back to a region of type 2.), either none or both the solutions which go from being complex 
conjugate to negative and back to complex conjugate must have acquired a min us si gn for their square root. This is 
due to the fact that in regions of type 2.), i.e. for d > 0, the right hand side of (A5) is always positive. To see this, 
consider the converse proposition, namely that AAP — 3iV < implies d < 0. First, note that 4M^ — 3A^ < ^ > 
and AlVP -3N > -3N. Check that at iM^ -3N = and at AM"^ -3N = -3N, d is manifestly negative; furthermore, 

3iV and finding the lone extremum at 4Af^ — 37V = ~N — 2N^/L, one verifies 



regarding ( A£ ) as a function of 4M^ 
that also there, d < 0, completing the argument. 

Apart from sign changes in acc orda nce with condition ( [A5| ), there is one other possibility for a continuous switch 
between the different solutions in ( A4), namely that two of these solutions may become degenerate (this is the generic 
case; higher degeneracies occur only at exceptional points in parameter space). This happens precisely when d = 0, i.e. 
when there is a concomita nt c hange of region. One can make this more precise: When starting with the conventional 
choice of square roots in (A4) in a region of type 1.), upon reaching a region of type 2.), Xi always remains unique, 
only two of the other three solutions may become degenerate. This remains true as long as one only alternates between 
regions of type 1.) and type 2.). On the other hand, if one continues from a region of type 2.) to a region of type 3.) 
(and back), bec omes degenerate with one of the other solutions; if, without loss of generality, one denotes as yi 
the solution of (A_2) which has remained real and positive, then it is X4 and 3:3 which become degenerate. This is the 
same switch which may also be effected by condition (A^) as shown further above. 

Armed with these properties, one can now argue that, despite all switches of branch which may happen as one 
tracks the solution through parameter space, it is only in regions of type 3.) that it may have a nonzero imaginary 
part (which up to trivial factors gives the sought-after eigenvalue distribution of the Hamiltonian (||)). As elucidated 
above, at ^ ^ 00 one starts in a region of type 1.) with the solution X4 in (A4). Now, changes of region occur precisely 
where d = 0. This latter condition is a fifth order equation in ^, so there can be at most five changes of region as ^ is 
varied, i.e. the following sequences are possible^: 

a. type (1) ^ type (2) -> type (1) -> type (2) -> type (1) ^ type (2) 

b. type (1) ^ type (2) ^ type (1) ^ type (2) ^ type (3) type (2) 

c. type (1) ^ type (2) ^ type (3) ^ type (2) ^ type (1) ^ type (2) 

d. type (1) ^ type (2) ^ type (3) ^ type (2) ^ type (3) ^ type (2) 

Of course, the sequences may be shorter if d has less than five real zeros as a function of ^. Now, according to the 
above, in sequences a.), b.), and d.) one always remains with solutions X4 and X3 (with the convention that yi is the 
solution which has remained real and positive throughout, otherwise 0:3 is replaced by one of the other two solutions). 
However, both in X4 and X3 in regions of type 2.) the imaginary parts of y/y2 and ^/ys cancel; thus one indeed only 
has a nonzero imaginary part of the solution in regions of type 3.). The remaining sequence c.) is more complicated. 
There, one may emerge in the second region of type 1.) with the solution 2:3 and thus, going to the last region of 
type 2.), switch to a solution in which the imaginary parts do not cancel anymore. Here, this last scenario will not 
be analyzed in detail. In practice, it suffices to note that the region of type 3.) already saturates the normalization 
of the eigenvalue distribution and that therefore there can not be any additional nonzero contribution from the last 
region of type 2.). 

Summing up, the imaginary part of Gho{0 is only nonzero in regions of type 3.), i.e. formally one can write 



1 



(All) 



where denotes the step function. The remaining ambiguity will be settled presently using condition (A_5) and the 
positivity of the eigenvalue distribution. To this end, one must specify more explicitly the solutions Wi of (A7). In 
the case c? < 0, which is the only one of interest here, (A7) is solved by 



Wi — 2Re (m) 



W2 



2Re (ue2"/3) W3 = 2Re (ue"^'^*/^) 



2iri/3\ 



(A12) 



^Note that a direct transition between regions of type 1.) and type 3.) can only occur if 60 ~ bi 
of ^, which is only possible at the exceptional point k = 0. 



in (A2) for some value 
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with 



Vd 



1/3 



(A13) 



Choosing the roots in ( [A13| ) in the conventional way (i.e. %/— T = i, any other choice merely permutes the Wi), one 
has Wi > W3 > W2, and consequently j/i > > 2/3 > 2/2, i-e- Ini^/yi = 0, the convention which was already used 
above. Furthermore, upon entering a region of type 3.) from a region of type 2.), an d ^/ y^ must be chosen to 
lie in different half-planes (upper or lower) in accordance with continuity and condition (A5), the right hand side 
of which is positive at the interface, AM^ — 3NM > 0. Only when within the region of type 3.), 4M^ — 3NM 
becomes negative, the root which touches zero acquires an additional minus sign; this is y/jfy according to the above 
conventions. Therefore, one can now specify IthGho to an overall sign. 



and this last ambiguity is resolved using the positivity of the eigenvalue distribution, which thus reads 
p//o(0=--ImGH„(0 



1 

( 

27r 



d)il - 9{-b2)e{bi)) v%^+ (1 - 20(4M3 _ 3A^Af))v^ 



(A14) 

(A15) 
(A16) 



APPENDIX B: EIGENVALUE DISTRIBUTION IN THE POLAR REPRESENTATION 

In the polar representation, expectation values 

(0> = ly d^rjJ{rf)0{rJ)eM~Nf3rf') 

= / d^r, 0(7?) exp I -iV/3 Y.'J^ +J2 " ''f )' + E | 



(Bl) 
(B2) 



are dominated by the saddle points of the exponent for large N. Note that the last term in the exponent in eq. 
(B2) is only of order 0{N) as opposed to 0{N'^) for the other two terms. Thus, going to a continuous eigenvalue 
distribution^ PBiv) ^01 N ^ 00, the saddle point distribution is determined by 







SpBiv) 



dxdy pB{x)pB{y)ln{x - y ) - P I dxxpsix) 



= J dxpB{x)Mx^-Tff-(3r]'^ 
or, taking the derivative with respect to rj and dividing by r/, 

= 2/^x4^+/? 



(B3) 
(B4) 

(B5) 



Here, the (as yet undetermined) edges of the support of ps have been made explicit. Note that c and h are positive 
because the eigenvalues rii are by definition positive. Substituting x"^ — s and if' = t, one obtains 



(B6) 



^The subscript makes explicit that pB is the eigenvalue distribution of the matrix B with eigenvalues rji. Below, also pg2 (x) 
ps(v^)/2%/S will be considered. 
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After defining 



Pb^{s) ^ PB^is)\l y^—^ (B7) 



one can use the inversion formulae on a finite interval |2C 



Jc^ y y ~ y ~ X TT J ^2 y h"^ — y X ~ y 



to obtain 



Pb^{s) = — pB2{s) = —d 2 PbW^ — \— 2 (^9) 



Now, the condition that ps^ must be normalized to one can be used to determine b, 

b^ = ^ + c' (BIO) 

On the other hand, the lower bound c is not determined by any auxiliary condition and should also be varied when 
determining the saddle point distribution. In other words, up to now, not the full space of possible ps has been 



explored, but only the subspace with fixed, albeit arbitrary, lower bound c. To determine c, one inserts the form (89) 



into the "action" in the square brackets in (P3) and now varies with respect to c. The first term in the action can 



be made manifestly independent of c by again substituting as in eq. (B6) and a shift in the integration variables; 
thus one very easily obtains c = 0, as one would expect. Finally, inserting the radius of the original semicircular 
distributions of the gauge fields, — 2/(3, one has 



PB^{s) = ^xl- (Bll) 

on the interval [0, 2a^] for the eigenvalue distribution of B^. 

APPENDIX C: SOLUTION OF FOURTH ORDER EQUATION DETERMINING EIGENVALUE 

DISTRIBUTION OF PARTICLE WITH SPIN 

The solutions of the equation (|8^ ) 

X* + C2a;^ + cix + Co = (CI) 

with 

C2 = -^{2L'^ + AKL + 6K'^ + 1) (C2) 
64 

ci = -^(2^^ + 12KL^ + ISK^L + 3L - 9K) (C3) 
16 

Co = - — {S2KL^ + ASK^L^ + 8L^ - A8KL + 9) (C4) 

are determined by the solutions of the corresponding cubic resolvent equation 

2/3 + 62y2 + biy + 6o = (C5) 

with 

16 

b2^- — il + 2K^ + L^ + {2K + Lf)<0 (C6) 
o 

2'i6 

6i = —(3 + + {3K - Lf + 24L^{K + L/Zf + {ZK + Lf/Z) > (C7) 
4096 

bo = -^^(2^^ + 12^i^ + ISA'^L + ZL- 9Kf < (C8) 
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as in eq. ( A4). Additionally, the signs of the square roots of the yi in ( A4) must be adjusted so as to guarantee 



VmVy^Vm = ci = -^(2i' + i2kl^ + isk^l + 3l - 9K) 



(C9) 



The third order equation (C5) in turn can be simplified by the substitution y = w — 62/3 to 



w + pw + q = 



(CIO) 



with 



256 

p= (2 - 3K^ - 9K^ - UKL - UK^L + + 2K^L^ + AKL^ - L") 

27 



q = 



4096 
729" 



(7 - 36K^ + 27K^ + 5AK^ + A8KL + lUK^L + lOSA'^i + 12L^ 



(Cll) 



(C12) 



and the associated discriminant 



27 ^ T 
4194304 

19683 



(3 - 24K^ + 54if4 - 81K^ + 8L^ + 260K^L^ 



2 r6 



+180if''i^ + - 128ii:^i* - 118if'*i^ + 20^^^^" - V) 



(C13) 



(CM) 



The analysis of these equations turns out to be considerably simpler than for the equations describing the case of the 
spinless particle. There is only one property which is not immediately obvious, namely that the discriminant (C14) 
has at most two zeros as a function of L for L G [0, 00] (remember that the whole problem is symmetrical about 
L = and it is therefore sufficient to consider positive L). One proves this by regarding d as a polynomial in 
and constructing the Sturm chain ||l^ of d{L'^) at = and for ^ 00. The absolute difference between the 
number of sign changes occuring in the one or the other chain gives the number of zeros of d{L^). One can check that 
this is less than or equal to two for any K, where it turns out to be advantageous to distinguish between the cases 
K"^ e [0,l/60],/i:2 g [1/60,1/3], and K"^ G [1/8,00]. 

Now, just as in the case of the spinless particle (cf. Appendix ^), the solutions of the cubic resolvent equation (C5) 
can in principle display three qualitatively different types of behaviour, partitioning the space of parameters K and L 
into regions corresponding to the different types. Here, however, there are no regions of type 3.), since one always has 
62 < and hi > 0. Therefore, the support of the eigenvalue distribution of the Pauli Hamiltonian must come from 
a region of type 2.) in which the imaginary parts of the two complex conjugate solutions do not cancel. Moreover, 
due to the fact that d has at most two zeros as L is varied over all positive values, and d < for L — > c», there is 
only one unique region of type 2.) for L > (and its mirror image for L < 0). Denoting without loss of generality as 
yi and 2/2 the two complex conjugate solutions and as j/3 the real positive solution, the form of the relevant solution 
X = aGQ^{0) — 4:L/3 of (|Cl|) on the region of type 2.) is therefore already almost uniquely determined: 



Ima; = --|^(^/y^ - y^) 



Rex = ±-y/y3 



(C15) 



where in liax it has additionally been used that due to the positivity of the eigenvalue distribution, Imx — aim Gq^ (0) 
must be chosen such that it is negative (cf. eq. (^0|)). The ambiguity in Rex on the other hand is resolved with the 
help of condition (|C9|). When the right hand side of ( |C9D is positive, it is consistent to choose the square roots ^/yi 
in the conventional way and ( C15 ) c an th en be identified with one of the two solutions X2 or 2:3 in (A4), i.e. 
then comes with the positive sign in (C15). On the other hand, if a sign change occurs in condition (C9) within the 
region of type 2.), then it can only be the real positive solution y^ which touches zero, i.e. it is which acquires 
an additional minus sign. In other words, one has 



Re x 



1{20{9K 



2L^ - 12KL^ - ISK^L - 3L) - 1)^^ 



(C16) 



It remains to give the solutions Wi of equation (CIC) explicitly. In a region of type 2.), i.e. for d > 0, they are 

W3 = u + u (C17) 



W2 - e-^^'/^u + e2-/3y/ 
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where 



IpI 



1+^ 



1/3 



M 

3 



-1/3 



(C18) 



This fixes the solutions yi — Wi — h^j'i of the cubic resolvent equation (|C5|), in terms of which the solution for Gq^ (0) 
is now finally given as 



ImGQ,(0) 



2a 



0(d)|z(Vyr-v^)i 



Re Gq^ (0) = ;^(20(9is: - - 12XL2 - ISisT^L - 3L) - 1)^^ + ^ for d > 



(C19) 
(C20) 
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FIG. 1. 



(A) for different momenta fc, in D = 3 dimensions. 



FIG. 2. p /7^(A) for different dimensions D, with momentum k — 0. 



FIG. 3. p /7i-(A) for different dimensions D, with momentum fc = 0, rescaled corresponding to keeping a fixed norm of the 



vector field as D is varied (see text). 



FIG. 4. Left-hand graph: Loci of solutions of det(QA — fi) — for a generic realisation of the gauge matrices Ai. Changing 
the momentum k, which enters the above equation as a parameter, merely amounts to shifting the origin of the A — /i-plane in 
the direction of the line A = —p/2 (see text). Right-hand graph: Non-monotonous behaviour of trajectories in the A — /i-plane 
is impossible (see text). Dotted lines serve merely to guide the eye. 



FIG. 5. p^/j^{X) for different momenta k. The reader is reminded that Hs was a 2A'' x 2A'^ matrix, and that that therefore 
the eigenvalue distribution of \/Hs as constructed in the text is normalised to two. 
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FIG. 6. P^j^ (solid lines) and (dashed lines) in three dimensions for (from left to right) k/a = 0, i, 1, 2, 5 (note the 

shift in the A-ax;is for the last two plots). 
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